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We analytically calculate the influence of a plasma on the shadow of a black hole (or of another 
compact object). We restrict to spherically symmetric and static situations, where the shadow is 
circular. The plasma is assumed to be non-magnetized and pressure-less. We derive the general 
formulas for a spherically symmetric plasma density on an unspecified spherically symmetric and 
static spacetime. Our main result is an analytical formula for the angular size of the shadow. As 
a plasma is a dispersive medium, the radius of the shadow depends on the photon frequency. The 
effect of the plasma is significant only in the radio regime. The formalism applies not only to black 
holes but also, e.g., to wormholes. As examples for the underlying spacetime model, we consider the 
Schwarzschild spacetime and the Ellis wormhole. In particular, we treat the case that the plasma 
is in radial free fall from infinity onto a Schwarzschild black hole. We find that for an observer far 
away from a Schwarzschild black hole the plasma has a decreasing effect on the size of the shadow. 

The perspectives of actually observing the influence of a plasma on the shadows of supermassive 
black holes are discussed. 

PACS numbers: 04.20.-q - 98.62.Sb - 98.62.Mw - 98.35.Jk 


I. INTRODUCTION 

In the last twenty years increasing evidence has been 
found for the existence of a supermassive black hole at 
the center of our galaxy. It is expected that a distant 
observer should “see” this black hole as a dark disk in 
the sky which is known as the “shadow”. It is some¬ 
times said that the shadow is an image of the event hori¬ 
zon, and the size of the shadow is estimated by the an¬ 
gle under which the observer would see the horizon ac¬ 
cording to Euclidean geometry. Actually, the boundary 
of the shadow corresponds to light rays that asymptot¬ 
ically approach the photon sphere (at r = 3 M in the 
Schwarzschild case) and not the horizon (at r = 2 M in 
the Schwarzschild case). Moreover, light rays do not fol¬ 
low straight lines in Euclidean geometry but they are 
bent. For these two reasons, the angular diameter of 
the shadow is actually bigger than the naive Euclidean 
estimate suggests. For the black hole at the center of 
our galaxy, it amounts to about 53 ^as whereas the Eu¬ 
clidean estimate gives only about 20 pas. At present, 
two projects are under way to observe this shadow which 
would give important information on the compact object 
at the center of our galaxy. These projects, which are go¬ 
ing to use (sub) millimeter VLBI observations with radio 
telescopes distributed over the Earth, are the Event Hori¬ 
zon Telescope (http://eventhorizontelescope.org) and the 
BlackHoleCam (http://blackholecam.org). 

On the theoretical side, the shadow is defined as the 
region of the observer’s sky that is left dark if there are 
light sources distributed everywhere but not between the 
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observer and the black hole. For a non-rotating black 
hole, the shadow is a circular disk. For a Schwarzschild 
black hole the angular diameter of the shadow was cal¬ 
culated, as a function of the mass of the black hole and 
of the radius coordinate where the observer is situated, 
by Synge 0- (Synge did not use the word “shadow”. 
He calculated what he called the “escape cones” of light. 
However, the complement of the opening angle of the es¬ 
cape cone at the observer gives exactly the diameter of 
the shadow.) For a rotating black hole the shadow is no 
longer circular. The first to correctly calculate the shape 
of the shadow of a Kerr black hole was Bardeen [2|. More 
generally, the shape and the size of the shadow can be 
calculated analytically for the whole class of Plebanski- 
Demianski spacetimes 0 0] • All these calculations are 
based on the assumption that light propagates along 
lightlike geodesics, without being influenced directly by 
matter. If such an influence is to be taken into account, 
one usually resorts to numerical calculation. In particu¬ 
lar, Falcke, Melia and Agol [g] have numerically simulated 
the visual appearance of the black hole at the center of 
our galaxy, assuming that it is a Kerr black hole, with 
scattering and the presence of emission regions between 
the observer and the black hole taken into account (at 
0.6 and 1.3 mm wavelengths). Sophisticated ray trac¬ 
ing programs have been written for producing realistic 
images of a black hole surrounded by an accretion disk, 
e.g. for the movie Interstellar. The numerical techniques 
used for this movie are described in detail, along with a 
review of earlier work, by James et al. 

While the purely geometric calculation of the shadow 
can be done analytically, at least for spacetimes in which 
the equation for lightlike geodesics is completely inte¬ 
grate, virtually all of the work that takes the influence 
of matter into account was based on numerics so far. In 
this paper we want to take the first steps towards an 
investigation of the shadow in matter based on analyti- 



2 


cal calculations. It is clear that this can be done only for 
highly idealized situations, but nonetheless we think that 
it is worthwile to do so. In contrast to numerical simu¬ 
lations, which can depict the situation only for a partic¬ 
ular choice of the parameters involved, analytical results 
demonstrate how exactly an effect depends on these pa¬ 
rameters and they bring out the general features more 
clearly. Also, analytical results are useful for testing nu¬ 
merical codes. In this paper we begin with the simplest 
non-trivial case: We want to analytically calculate the in¬ 
fluence of a non-magnetized pressure-less plasma on the 
size of the shadow of a non-rotating black hole. 

A plasma is a dispersive medium, i.e., the light rays 
deviate from lightlike geodesics in a way that depends on 
the frequency. The Hamiltonian for the light rays can 
be derived from Maxwell’s equations where the source of 
the electromagnetic field consists of two charged fluids, 
one modelling the ions and the other the electrons. For 
a plasma on a curved background, the transition from 
Maxwell’s equations to ray optics has to be done by a 
two-scale method. A rigorous derivation of the Hamil¬ 
tonian for light rays was provided by Breuer and Ehlcrs 
Ij7i, Jj who considered a magnetized pressure-less plasma. 
For the much simpler case of a non-magnetized pressure¬ 
less plasma, a similar derivation can be found in Per- 
lick [9|. In the latter case the equation of light rays 
can be fully characterized by a scalar, i.e. direction- 
independent, index of refraction which is a function of 
the spacetime point and of the frequency. Then the re¬ 
sulting equation of motion for light rays belongs to a 
class that was studied in the text-book by Synge 0- 
The resulting formula for the deflection angle of light in 
a plasma whose density is a function of the radius co¬ 
ordinate was calculated in the Schwarzschild spacetime 
(and more generally in the equatorial plane of the Kerr 
metric) by Perlick [§]. The same result was found and 
further discussed by Bisnovatyi-Kogan and Tsupko 0, 
0,0. Morozova et al. 0 generalized the calculation 
to the Kerr metric off the equatorial plane, assuming that 
the black hole is slowly rotating. For recent suggestions 
of how the effects of a plasma on light rays could actually 
be observed we refer to Er and Mao jT^ and to Rogers 
0 . 

In this paper we calculate the angular diameter of the 
shadow for a non-rotating black hole that is surrounded 
by a non-magnetized pressure-less plasma. Although we 
are mainly interested in the Schwarzschild case, in large 
parts of the paper we work in an unspecified spherically 
symmetric and static spacetime, so the results can be 
applied, e.g., also to charged black holes or to wormholes. 

In our work the plasma manifests itself as a refractive 
and dispersive medium, which leads to a change of the 
light rays in comparison with vacuum. As a result, the 
geometrical size of the shadow is changed. We do not take 
into account the processes of absorption and scattering 
of the photons by plasma electrons. Also, we neglect the 
gravitational field of the plasma particles. 

The paper is organized as follows. In Section [TT] we 


work out the equation of motion for light rays in a plasma 
on a spherically symmetric and static spacetime. In Sec- 
tion lllll we determine the circular light orbits which are of 
crucial relevance for the formation of the shadow. In Sec¬ 
tion [TV] we calculate the angular diameter of the shadow. 
We specify the results in Section [V] to the case of a 
low-density plasma and in Section IVII to the case that 
the plasma has the particular density profile that corre¬ 
sponds to spherically symmetric accretion of a dust onto 
a Schwarzschild black hole. In the same section we also 
discuss the perspectives of actually observing the influ¬ 
ence of a plasma on the shadow of a supermassive black 
hole. 

We use units such that G = c = 1, so the Schwarzschild 
radius is Rs = 2 M. Our choice of signature is 
{—,+,+,+}. We use Einstein’s summation convention 
for latin indices which take the values i 7 k, ... = t,r, d, <p. 


II. EQUATIONS OF MOTION FOR LIGHT 
RAYS IN A NON-MAGNETIZED PLASMA 


We consider a spherically symmetric and static metric 

gikdx l dx k = — A{r)dt 2 -\-B{r)dr 2 +D{r)( y d’d 2 +s\TL‘ n ddip 2 ^ , 

(!) 

where A(r), B(r) and D(r) are positive. We are mainly 
interested in the case of a Schwarzschild black hole but 
for the time being there is no need for specifying the 
metric any further. 

We assume that the spacetime is filled with a non- 
magnetized cold plasma whose electron plasma frequency 
uj p is a function of the radius coordinate only, 

A-irp^ 

Up( r ) = - N(r). (2) 

m 

Here e is the charge of the electron, m is the electron 
mass, and N(r) is the number density of the electrons 
in the plasma. The refraction index n of this plasma 
depends on the radius coordinate r and on the frequency 
uj of the photon as it is measured by a static observer, 


n(r, w) 2 = 1 — 



( 3 ) 


Because of the spherical symmetry we may restrict to the 
equatorial plane d = 7 t/2, p& = 0. Then the Hamiltonian 
for light rays in the plasma reads 

H = \(g lk PiPk+u p (r) 2 ) = (4) 


= I ( Pt + P r + P V 
2 ^ A{r) B(r) D(r ) 



For a derivation of the Hamiltonian (ED from Maxwell’s 
equations with a two-fluid source we refer to the litera¬ 
ture given in the introduction. 
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The light rays are the solutions to Hamilton’s equa¬ 
tions 


dH 


dH 


Pi = > x = o - 

OX 1 opi 


which in this case read 


dH 


P V = 


Pr = 


dH 

dp 

dH 

dr 


= 0. 


( 5 ) 


( 6 ) 


( 7 ) 


( 8 ) 


p 2 A'(r ) p 2 B'(r) p^,D'(r) d 


v> _ 

2 V A{r) 2 ' B{r) 2 ' D(r) 2 dr 


+ : 


oi r ) 2 ) , 


t = 


dH 

dp t 


Pt 


<P = 


dH 

dpu, 


dH 


A(r) 


P v 


D(r )’ 


Pr 


dp r B(r) ’ 


( 9 ) 


( 10 ) 


(ID 


with H = 0, i.e. 


0 = - 


Pt 


+ 


A(r) B(r) D(r) 


+ w p (r) 2 . (12) 


Here a dot means differentiation with respect to an affine 
parameter A, and a prime means differentiation with re¬ 
spect to r. 

^From (ED and 0 it follows that p t and p v are con¬ 
stants of motion. We write ujq '■= —pt- If w>o has been 
fixed and if A(r) —»■ 1 for r —> oo, which is the case if the 
spacetime is asymptotically flat, the frequency ui mea¬ 
sured by a static observer becomes a function of r by the 
gravitational redshift formula, 


w(r) = 


w 0 


VMr 


(13) 


By m, a light ray with constant of motion wo is re¬ 
stricted to the region where 


u 0 

A (r) 


> ixpir ) 2 


(14) 


The restriction (ED physically means that the photon 
frequency at a given point, w(r), must be bigger than 
the plasma frequency, w p (r), at the same point. This is 
always true for light propagation in a plasma. 


To derive the orbit equation, we use ED and ED to 
find 


dr 

dp 


r 

<P 


D(r)p r 


B{r)pp ' 

Upon substituting for p r from CC3b this results in 


dr_ _ ± \JD(r) 
dp y/B{r 



Zjh(r) 2 - 1 


where we have defined the function 


h(r) 2 = 


B{ r ) 

A{r) 


1 — A{r) 


u P (rf 


(15) 


(16) 


(17) 


In general, the orbit has to be decomposed into sections 
where r is increasing as a function of p and sections where 
it is decreasing, and the sign in (fTBl) has to be chosen 
appropriately. For a light ray that comes in from infinity, 
reaches a minimum at a radius R, and goes out to infinity 
again, integration over the orbit gives the formula for the 
bending angle, <5, 


■6 = 2 


r VW) 

Ir v d (. t ) 


^h(r) 2 - 1 

Ptp 


- 1/2 


dr. (18) 


As R corresponds to the turning point of the trajectory, 
the condition dr/dp\ R = 0 has to hold. This equation 
relates R to the constant of motion p v /u) o, 


h(R) 2 = . 


w 


(19) 


Then the deflection angle can be rewritten as a function 
of only R and Wq (for a given plasma distribution) as 


6 = 2 


VW) ( Hr) 


'r \H r )‘ 


- 1 


- 1/2 


dr. (20) 


Example 1: Schwarzschild spacetime 
For the Schwarzschild spacetime, 


A(r) = B(r) 1 = 1 — —— , D(r) = r 2 
r 


the function h{r) specifies to 
h(r) 2 = r 


2 _ 2 I T UpHf 


r — 2 M uj q 


Then the bending angle reads 


6 = 


( 21 ) 


( 22 ) 


( 23 ) 



, —!/2 


— 1 


dr 


yfr\/ r — 2 M 
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This formula for the bending angle in a plasma on 
Schwarzschild spacetime was derived in [9]. In [r| it was 
rederived using Synge’s approach and rewritten in terms 
of an elliptic integral for a homogeneous plasma; there 
also the strong deflection limit (6 O 1) was investigated. 

Example 2: Ellis wormhole 

As a second example we choose the Ellis wormhole §52} 
which is a traversible wormhole of the Morris-Thorne 
class 0. It is true that the existence of such wormholes 
is questionable because they need exotic matter 0] and, 
at least for a certain kind of perturbations, it has been 
shown that the Ellis wormhole is unstable 0-0|. On 
the other hand, wormholes have met with great inter¬ 
est because they make time travel possible and the Ellis 
wormhole is an instructive example for illustrating the 
applicability of our results. In this case the metric coef¬ 
ficients are 

A(r) = B(r ) = 1 , D(r) = r 2 + a 2 , ( 24 ) 

where the coordinate r ranges from — oo to oo and a is a 
constant that determines the radius of the “neck” of the 
wormhole. Note that this is an example of a spherically 
symmetric and static spacetime where we cannot make 
a transformation of the radius coordinate, r —> f, such 
that D(r) = f 2 . The reason is that the function D(r ) 
has vanishing derivative at r = 0 (i.e., at the “neck”), so 
the desired transformation fails to be a good coordinate 
transformation on any radius interval that contains the 
point r = 0. For the Ellis wormhole the function h{r) 
reads 

M0 2 = (r 2 + « 2 )(l-^) , (25) 

and the bending angle is given by 

7r + 6 = (26) 

9 f°° ( (r 2 +a 2 )(^g -u p {r) 2 ) _ ^ 1 /2 dr 
Jr \(R 2 +a 2 )(u, 2 ~u p (R) 2 ) ) 

III. CIRCULAR LIGHT ORBITS 

We now derive the condition for cicular light orbits 
which will be crucial for determining the shadow. Along 
a circular light orbit we must have f = 0 and r = 0. The 
first condition, by m, implies p r = 0; from m we get 
the equation 


0 = - 


w 


P% 


A[r) D(r) 


-Wp(r) 5 


On the other hand, ED implies 


(27) 


^From this equation we read that r = 0 together with 
r = 0 leads to p r = 0, and from (j8j) we get the second 
equation for circular light orbits, 


wgA'(r) P 2 p D'(r ) d 


° A(r) 2 ' + D(r ) 2 dr ^^‘ 


We solve each of these two equations m and 


(29) 

for 


P% = 




(30) 


P 


2 

V> 


D(r) 2 ( eo 2 A\r) 
D'(r) V A(r) 2 



(31) 


Subtracting these two equations from each other yields, 
after some elementary re-arrangements, the equation for 
the radius of a circular light orbit in the form 

0 = ±h(rf (32) 

with the function h{r) 2 from (1171) . Any solution r = r p h 
of (|32D determines the radius of a photon sphere. If a light 
ray starts tangentially to such a sphere it will stay on a 
circular path with radius r p h forever. If the spacetime is 
asymptotically flat, and if ui p (r) —> 0 for r —> oo, the out¬ 
ermost photon sphere is always unstable with respect to 
radial perturbations. This means that the circular pho¬ 
ton orbits in this photon sphere can serve as limit curves 
for light rays that approach them asymptotically. The 
radius r p h of the outermost photon sphere is the criti¬ 
cal value of the minimal radius R mentioned above. If 
a light ray comes in from infinity and reaches a mini¬ 
mum radius R bigger than r p h, it will go out to infinity 
again. The case R = r p h corresponds to a light ray that 
spirals asymptotically towards a circular photon orbit in 
the sphere of radius r p h- All other rays cross the photon 
sphere and we exclude the case that they can come back. 
The latter case can occur only if there is a second photon 
sphere. 

In the vacuum case, ui p (r) = 0, the condition for a cir¬ 
cular light orbit in an unspecified spherically symmetric 
and static spacetime was first given by Atkinson 0]. It 
is easy to check that our condition (l32l) reduces, indeed, 
to Atkinson’s if ui p (r) = 0. 


Example 1: Schwarzschild spacetime 

For the Schwarzschild spacetime, where h{r) is given by 

(H2), condition (1521) for a photon sphere reads 


r(r-3M) _ u p (r) 2 _ u p {r)u' p (r) 
(r — 2 M) 2 oig Wq 


(33) 


Pr = ~rr (B{r) r) = rB(r) + f 2 B'(r). 
dX 


For the special case that ui p depends on r via a power 
(23) law, this condition was already derived by Rogers [r|- If 
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there is no plasma, to p jr ) = 0, (|33l) gives the well-known 
result r p h = 3 M. For the sake of curiosity, notice that 
(13311 is identically satisfied if 


hjr ) 2 = r 2 


r 

r — 2M 


u P (r) 2 \ 
w o ) 


= C, 


(34) 


i.e. ijJ p (r) 2 


to n 


r 

r — 2M 



(35) 


with a positive constant C. For this particular density 
profile of the plasma, there is a circular photon orbit with 
constant of motion too at any radius r for which 


r 


3 


r-2M 


> C. 


(36) 


Example 2: Ellis wormhole 

For the Ellis wormhole, hjr) is given by (HT51) and condi¬ 
tion (l32l) for a photon sphere reads 


0 = 



u P ( r ) 2 \ 

w o ) 


jr 2 + a 2 ) 


U} p( r ) UJ P ( r ) 


(37) 


Without a plasma, to p (r) = 0, there is a unique photon 
sphere at the neck of the wormhole, r p h = 0. The same is 
true for a homogeneous plasma, to p (r) = constant. How¬ 
ever, in an inhomogeneous plasma there may be arbitrar¬ 
ily many photon spheres. For an observer at a large posi¬ 
tive r coordinate, the outermost photon sphere is relevant 
for the formation of the shadow which will be discussed 
in the next section. 


IV. RADIUS OF THE SHADOW 

The shadow of a spherically symmetric and static black 
hole is defined in the following way, see Fig. [T] Consider 
light rays sent from an observer at radius coordinate rp 
into the past. As we want to take the influence of a 
plasma into account, by a light ray we mean a solution 
to the equations of motion discussed in Section [Tl] These 
light rays can be divided into two classes: Light rays of 
the first class go to infinity after being deflected by the 
black hole. Light rays of the second class go towards the 
horizon of the black hole. If we assume that there are 
no light sources between the observer and the black hole, 
initial directions of the second class correspond to dark¬ 
ness on the observer’s sky. This dark circular disk on the 
observer’s sky is called the shadow of the black hole. The 
boundary of the shadow is determined by the initial di¬ 
rections of light rays that asymptotically spiral towards 
the outermost photon sphere. Here it is crucial that the 
light rays in the photon sphere are unstable with respect 
to radial perturbations because otherwise they could not 
serve as limit curves. For this reason, the construction 
of the shadow works for any spherically symmetric and 


static spacetime that admits an unstable photon sphere. 
This includes not only black holes but also, e.g., worm- 
holes. We have already mentioned that in a spherically 
symmetric and static spacetime that is asymptotically 
flat the outermost photon sphere is always unstable, pro¬ 
vided that the plasma density tends to zero for r — > oo. 
We will now calculate, for this situation, the angular ra¬ 
dius a s h of the shadow. The observer is assumed to be 
static somewhere between the outermost photon sphere 
and infinity. 

We consider a light ray that is sent from the observer’s 
position at ro into the past under an angle a with respect 
to the radial direction. ^From Fig. [T] we read that a is 
given by 


cot a = 


\/ fjrr dr 


y/glpip dip 


r=r o 


\jB(r) dr 
\JD(r) dip 


r=ro 


(38) 


If the light ray goes out again after reaching a minimum 
radius R , the orbit equation (1161) can be rewritten, with 
the help of (TTUl) . as 


dr_ = y/D(r) / h 2 jr) 
dip \/B(r) V h 2 (R) 

For the angle a we obtain 


Using 


cot 2 a = 


h 2 jrp) 
h 2 (R ) 


- 1. 


1 + cot 2 a = 


1 

sin 2 a 


we get 


sin 2 a 


h(R) 2 
h{r Q ) 2 ' 


(39) 


(40) 


(41) 


(42) 


The boundary of the shadow a s h is determined by light 
rays that spiral asymptotically towards a circular light 
orbit at radius r p h- Therefore the angular radius of the 
shadow is given by sending R — »• r p h in (l42l) , 


sin 2 a sh 


h{r ph ) 2 

h{r Q ) 2 ’ 


(43) 


where h(r ) is given by the formula (1171) . 

For many applications we may assume that the ob¬ 
server is in a region where the plasma density is negligibly 
small. Then m implies 


h(r 0 ) 2 


Dj.ro) 

^(t-q) 


and (l43l) reduces to 


( 44 ) 


sin 2 a sh 


Djr ph )Ajr 0 ) / 
Ajr ph )D(r 0 ) V 


Ajr p h)to p {r p h ) 2 \ 
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FIG. 1. For R —> r p h, the angle a approaches the angular 
radius a s h of the shadow. 


N(r 0 ) < N(r ph ). 

This demonstrates that, under the assumptions stated, 
the plasma always has a decreasing effect on the size of 
the shadow. 

We emphasize that the preceding calculation applies 
not only to black holes but also to other spherically sym¬ 
metric and static spacetimes with an unstable photon 
sphere, e.g. to ultracompact stars and to wormholes. As 
long as there is no light coming from the direction of 
the central object to the observer (which means, in par¬ 
ticular, that the central object must not have a bright 
surface), these objects would cast a shadow in the same 
way as a black hole. 

Let us summarize the results of this section. To find 
a s h for a given metric of the form m, a given plasma 
concentration iV(r), a given photon frequency at infinity 
uj o and a given observer position ro , we have to calculate 
r p h using ea. (l351) and to substitute the result into formula 
(®. Note that the photon frequency at the observer 
position is w(r-o) according to (1T31) . 


Example 1: Schwarzschild spacetime 
For the Schwarzschild case (H51) specifies to 


sm a sh = 


2 ( r ph 
?ph \ r p h — 2M 


.2 ( r Q 
° \r 0 - 2 M 


w P (r p h ) 5 


o(roY 


(46) 


where r p h has to be determined from (1531) . In Section lYll 
below we will evaluate this equation for a particular 
plasma density profile. For vacuum, w p (r) = 0, our con¬ 
sideration gives 


Krf 


1 - 2 M/r’ 


r ph = 3 M, 


. 2 27M 2 (1 - 2M/ro) 

snrdsh = -- —2 ---- • (47) 

r o 

This is Synge’s [l| formula for the radius of the shadow of 
a Schwarzschild black hole which was mentioned already 
in the introduction. 


Example 2: Ellis wormhole 

We now calculate the radius of the shadow for the Ellis 
wormhole. In this case, our assumption of no light com¬ 
ing towards the observer from the direction of the cen¬ 
tral object means, in particular, that there are no light 
sources in the region r < 0. For the Ellis wormhole, the 
function h[r) is given by (|25l) . so the equation (l43l) for 
the angular radius of the shadow specifies to 


_ ( r ph + Q2 ) P h) 2 ) 

( r O + ° 2 ) (^0 - Wp(r-o) 2 ) ' 


(48) 


For a homogeneous plasma, co p (r) = constant, we have 
r p h = 0 and thus 


sin 2 a sh = o ~— 5 -, (49) 

r o + 0,2 

so a homogeneous plasma has no influence on the size of 
the shadow. 


V. THE SHADOW IN A LOW-DENSITY 
PLASMA 


If the plasma frequency is much smaller than the pho¬ 
ton frequency, the equations for the photon sphere and 
for the radius of the shadow can be linearized about the 
corresponding values for vacuum light rays. To work this 
out, we rewrite (El as 


h(r) 2 = Q(r)(l - e/3(r)^ (50) 


where 


_ D(r) A(r)ujp(r) 2 uj p (r) 2 

A(r) ’ Wg w(r ) 2 


(51) 


Here we have introduced a book-keeping parameter e 
which will be set equal to unity after all equations have 
been linearized with respect to e. 

The equation for a photon sphere, (1351) . now reads 


0 = Q\r) (l - e/3(r)) - eQ(r)/3'(r). (52) 

We write the solution to this equation in the form 

e p h = 7" p h + er ph + ■ • ■ (53) 

where r|( h is a solution in the case without plasma, i.e., 

Q'( r ph) = 0 ■ (54) 
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After inserting these expressions into (1521) and comparing 
coefficients of e we find 


= 


Q"K h ) 


(55) 


Depending on the signs of /3'(r p h) and Q"(rp h ), this ex¬ 
pression can be positive or negative, i.e., the plasma can 
shift the photon sphere inwards or outwards. 

We now insert the expansion (1531) into the equation for 
the shadow (H51) . If we neglect all terms of quadratic and 
higher order in e, and set e equal to unity again, we find 


0(r° ) / \ 

sin 2 a sh = Q( r P p) ( 1 -/ 3 ( r ph)+/3(^o)) , w P {r) < w(r). 

(56) 

Note that Cp h has dropped out from the equation for a s h 
to within this order. According to dSHl) . the plasma has 
a decreasing effect on the shadow as long as /3(ro) < 

/^ph)- 

Let us summarize the results of this section. To find 
the radius of the photon sphere r p h for a given metric 
ED, a given photon frequency at infinity uj o, and a given 
plasma frequency that satisfies the condition w p (r) -C 
u(r) one has to proceed in the following way: Write down 
the functions Q(r) and /3(r), see (I5l1) : calculate r-p h from 
(1511) and Tp h from (l55l) : then calculate r p h according to 
(15H1) with e = 1. To find the angular radius a s h of the 
shadow for a given observer position ro in this case, we 
have to insert r° h into (l56l) . 


Example 1: Schwarzschild spacetime 
As an example of the approximation formalism developed 
in this section, we consider the Schwarzschild spacetime 
for the case that the plasma electron density is given by 
a power law, 




(57) 


where /3q and k are positive dimensionless constants. 
Then r° h = 3 M and 


i _ /3o M 

r ph — gfc+l 



(58) 


i.e., depending on the density profile the radius of the 
photon sphere may become smaller (k > 2) or bigger 
(k < 2) than 3 M. In the case that k = 2 we obtain 
r P h = 0- Note that in this case the equation for the 
photon sphere can be solved exactly and gives the same 
result, r p h = 3 M. The first-order equation for the radius 
of the shadow, (15H1) , yields 

sin 2 a p h = (59) 



If the observer is far away from the black hole, ro M, 
this can be simplified to 


. 2 27 M 2 ( 

sin a p h = — 2 — 1 

r o V 



(60) 


Example 2: Ellis wormhole 
For the Ellis wormhole, 

Q(r)= r 2 +a 2 , /3(r) = . (61) 

With Tp h = 0, the first-order approximation yields 

i _ a 2 M°K(°) 


sin a sh = 


( r o + « 2 ) 


a; p (0) 2 M rp) 2 \ 

w o w o ) 


(62) 


VI. SPHERICALLY SYMMETRIC ACCRETION 
OF A PLASMA ONTO A SCHWARZSCHILD 
BLACK HOLE 

In this section we consider in greater detail the special 
case that the underlying spacetime is the Schwarzschild 
spacetime and that the plasma electron density corre¬ 
sponds to spherically symmetric accretion of dust. As 
usual, we use the word “dust” as a synonym for a 
pressure-less perfect fluid. Vanishing pressure implies 
that the flow lines are geodesics, i.e., that the fluid parti¬ 
cles are freely falling. Note that for our plasma model the 
Hamiltonian and, thus, the equations of motion for light 
rays do not depend on the velocity of the plasma. This is, 
of course, different for other types of media, see Synge’s 
book E3. In the case of an infalling plasma we have to 
determine the plasma density which is the only quantity 
that enters into the equation for the plasma frequency. 

We start out from the continuity equation 


9i(V-gpu l ) = 0. (63) 

Here p is the rest-mass density of the plasma, u l is the 
4-velocity, and g is the determinant of the metric. For 
the sake of simplicity, we consider a neutral hydrogen 
plasma and we assume that the electrons have the same 
4-velocity as the protons. As the electron mass is negli¬ 
gibly small in comparison to the proton mass, p can be 
written as p = m p N where m p is the proton rest mass 
and N is the number density of the protons (which, as 
the plasma is assumed to be neutral, coincides with the 
number density of the electrons). 
















For spherically symmetric and stationary accretion in 
the Schwarzschild spacetime the continuity equation re¬ 
duces to 

— (r 2 p(r)u r (r)) = 0 . (64) 

dr 

Integration gives 

4nr 2 p(r)u r (r) = —Ma = const (65) 

where Ma is a stationary mass flux. 

We assume that the particles are dropped in radial free 
fall from rest at infinity, i.e., that our plasma is so “cold” 
that the pressure can be neglected. Then the integral 
curves of u M are radial geodesics, so the radial component 
of the 4-velocity is (2lj 


dr [2M 
~dr~~\ 


r 


( 66 ) 


Thereupon, eq. (Iti5l) gives us the rest-mass density, 


sina,i,(f3 0 ) 



FIG. 2. Dependence of the radius of the shadow on /3o for 
different observer positions ro: ro/M = 3.333 (a), 3.768 (b), 
5 (c), 10 (d), 50 (e). The dependence is calculated using the 
exact formulae. It is shown that for a distant observer the 
radius of the shadow becomes smaller. 


^ M a 1 
P[l) 4 ;rVlMrW 


(67) 


Equations (l66l) and (l67l) also follow from Michel’s pio¬ 
neering paper [21 ] on fully relativistic spherically sym¬ 
metric accretion where we have to specify Michel’s equa¬ 
tions (9) and (10) to the case of a dust. 

The rest-mass density p(r ) gives us the ratio of fre¬ 
quencies 


u} p (r) 2 4ne 2 N(r) 47re 2 p(r) 

Wq rn e ojQ m e m p ujQ 


M 3 / 2 

r 3/2 


( 68 ) 


As an aside, we mention that the 4-velocity component 
u r directly gives the velocity of the infalling particles as 
measured by a static observer. To demonstrate this, fol¬ 
lowing [2Cj (see also @|), we first observe that (l66l) to¬ 
gether with the normalization condition giju l u° = — 1 
implies that 


hence 



(71) 


where 


2 M a 


m e m p ujQ\/2Al M 3 / 2 


(69) 


For this function w p (r), the angular radius (1461) of the 
shadow is plotted, as a function of /3o for different ob¬ 
server positions, in Fig. [2] 

Assuming that the approximation of a low-density 
plasma is justified, we can use formula (1591) and find 


sin 2 a sh = 


(70) 


dr dr dr 12M ( 2 M\ 

dt dr dt V r \ r ) 


(72) 


For a static observer at radius coordinate r, a coordinate 
time interval dt corresponds to a proper time interval 


dt' = 



(73) 


On the other hand, a radial coordinate separation dr cor¬ 
responds to a proper radial distance equal to 


27 M 2 / 

1 2M ) 


2 M\ 

iA)M 3/2 \ 

dr' = ( 

i-“) 

r 2 V 

r o v 

ro ) 

ro J 

1 4 /2 ) ' 

V 

r J 


dr. 


(74) 


In this case the function — /?(rp h ) + /3(ro) which enters 
into is equal to zero when ro — 3.768 M and it takes 
its maximum when ro = 3.333 M. When ro > 3.768M, 
the function is negative. Therefore, if the approximation 
of a low-density plasma is justified, the shadow becomes 
bigger if ro < 3.768 M and it becomes smaller if ro > 
3.768 M. 


As a consequence, the velocity of the radially infalling 
particle, as measured by a static observer at r, is given 
by 



(75) 


We will now use the results of this section for estimating 
the numerial values of the plasma effects on the shadow of 
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the supermassive black holes at the center of our galaxy 
and of M87. 

For these applications we can always use the formulas 
with ro Rs, so the angular radius of the shadow (IfUl 
can be written as 

. 2 27 Rg A) tr! o\ 

sh = (1 “ Z = 3572 ■ ( 70 ) 

/3o was defined in (jfifljl . For estimating the plasma correc¬ 
tion numerically, we re-express all quantities in Gaussian 
cgs units, i.e., we restore factors of c and G. Then /3 0 
reads 


e 2 MAC 3 


m e m p uj q v / 2GM(GM) 3 /2 ' 


(77) 


The mass accretion rate Ma can be estimated via the 
observed luminosity of the galactic center. We write 

L ~ r,M A c 2 , (78) 


where 77 is a non-dimensional coefficient characterizing 
the accretion efficiency. We obtain for the correction due 
to the presence of the plasma: 


2 3 / 2 e 2 L 1 
3 5 / 2 m e m p c rjc 2 


(79) 


Note that the circular frequency ui 0 is connected with the 
ordinary frequency uq by the usual relation, wo = 2ttvq 1 
and that in the expression for Z we can assume, as a valid 
approximation, that u>o is related to the wavelength Ao 
by the vacuum relation w 0 = 27tc/Ao. 

The value of 77 depends on the specific model of accre¬ 
tion. For estimations we consider 771 = 10 -4 (see 0 . 
0) and 772 = 0.1 (see 0]). 

Let us consider the galactic center source Sgr A*. For 
estimations we put M = 4.3 x 1O 6 M 0 , ro = 8.3 kpc (see 
0 , 0 ). Then Synge’s formula (l47ll gives for the angu¬ 
lar radius of the shadow in vacuum a value of a s h ~ 27 
7 ms, i.e., a diameter of about 53 ^tas which is expected 
to be resolvable with Very Long Baseline Interferometry 
(VLBI) soon 0, 0. By contrast, a naive Euclidean 
estimate of the horizon size, a ~ tana = Rs/to, yields 
a ~ 10 //as, i.e., a diameter of about 20 //as, as was 
already mentioned in the introduction. For a rotating 
(Kerr) black hole, the shadow is flattened on one side, 
but its vertical diameter is still given by Synge’s formula 
if ro » Rs, see 0. 

For estimating the correction Z we use L = 10 6 L 0 
(see, for example 0 ). For Ao = 1 mm we obtain: 

Z = 0.8 x 10 -5 for 77 1; Z = 0.8 x 1CT 8 for r) 2 . (80) 


Note that observations of the shadow are planned in the 
(sub)millimeter regime because at wavelengths of more 
than about 1.3 mm it is expected that the shadow is 
washed out by scattering @. We see that in this regime 
the effect of a plasma is rather small, for the chosen values 


of 77 . However, for bigger radio wavelengths the effect can 
be significant. For Ao = 10 cm we obtain: 

Z = 0.8 x 1CT 1 for 771 , Z = 0.8 x 10 " 4 for 772 . (81) 

We now turn to the galactic center source M87. For 
estimations we put M = 3 x 1O 9 M 0 , ro = 18 Mpc, 
L = 7 x 10 4 ° ergs s _1 (see, for example 0). Then 
Synge’s formula leads to an angular radius of the shadow 
in vacuum of a s h — 9 //as. 

For Ao = 1 mm we obtain: 

Z = 0.3 x 10 -9 for 771 , Z = 0.3 x 10 " 12 for 773 . (82) 

For Ao = 100 cm we obtain: 

Z = 0.3 x 10 " 3 for 771 , Z = 0.3 x 10 -6 for 772 . (83) 


VII. CONCLUSIONS 

In this paper the first steps towards an investigation 
of the shadow in the presence of matter based on analyt¬ 
ical calculations have been taken. We have analytically 
calculated the influence of a non-magnetized pressure¬ 
less plasma on the angular size of the shadow. It was 
our goal to derive all relevant formulas for an unspeci¬ 
fied spherically symmetric and static spacetime, before 
specializing to the case of a Schwarzschild black hole. 
The gravitational field was not supposed to be weak. We 
have worked in the frame-work of geometrical optics; in 
this approximation the presence of the plasma leads to a 
change of the geometrical size of the shadow via a change 
of the light rays in this medium. 

The equation of motion for light rays in a plasma with 
a spherically symmetric density distribution on a spher¬ 
ically symmetric and static spacetime © has been de¬ 
rived, see (1161) . In particular, the formula for the photon 
deflection angle has been shown for unbound orbits, see 
We have also found the radii of the circular light 
orbits (radii of the photon spheres), see (13^1) . The cen¬ 
tral result of our paper is an analytical formula for the 
angular size of the shadow, see m- □ 


1 Soon after a preprint of the present paper was made public on 
the arXiv flTlL a paper by Atamurotov and Ahmedov fl8ll ap¬ 
peared, again on the arXiv, where the shadow was calculated for 
the case of a homogeneous plasma. Their results are in obvious 
disagreement with ours if the latter are specialized to the case of 
a homogeneous plasma. The reason is in the fact that they erro¬ 
neously treated the index of refraction as a constant. Actually, 
the index of refraction involves the plasma density, which is con¬ 
stant in a homogeneous plasma, and the photon frequency, which 
in a gravitational field depends on the space point. Therefore, the 
index of refraction is not a constant for a homogeneous plasma. 
This has the consequence that in the Schwarzschild spacetime 
the photon sphere in a homogeneous plasma is different from the 
vacuum case. In the meantime, a strongly revised version of their 
paper appeared in PRD fl9ll where they are now considering a 
plasma density proportional to r -1 . 







10 


Special attention was given to the realistic case when 
the plasma frequency is much smaller than the photon 
frequency, see Section El We have shown that in this 
case the plasma has a decreasing effect on the size of 
the shadow provided that the plasma density is higher 
at the photon sphere than at the observer position, see 
(1561) . In particular, a simple formula for the size of the 
shadow was presented for the case of a power-law density 
distribution and an observer at a large distance, see ffl- 
In the presence of a plasma the size of the shadow 
depends on the wavelength at which the observation is 
made, in contrast to the vacuum case where it is the 
same for all wavelengths. For the underlying space- 
time we have treated two examples, the Schwarzschild 
black hole and the Ellis wormhole. In particular, the 
case of spherically symmetric accretion of plasma onto 
a Schwarzschild black hole was considered in detail, see 
Section IVII We have found that for an observer far away 
from the Schwarzschild black hole the plasma makes the 
shadow smaller. As examples, we have considered Sgr A* 
and M87. For the specific accretion model used here we 
have found that the effect of the presence of a plasma on 
the size of the shadow can be significant only for wave¬ 
lengths of at least a few centimeters. At such wavelengths 


the observation of the shadow is made difficult because 
of scattering. 

Apart from considering more complicated plasma mod¬ 
els, an obvious next step would be to generalize the anal¬ 
ysis to the case of axially symmetric and stationary situ¬ 
ations to include rotating black holes. Then the shadow 
is no longer circular, and the plasma has an effect not 
only on its size but also on its shape. We are planning 
to work out the details in future work. 
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